On the origin of the rise of F2 at small x 
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Abstract 

We show that, provided that the non-perturbative input is regular at the 
right of the uj = singularity of the dominant DGLAP anomalous dimension, 
the rise of F2 at small x, experimentally measured by the averaged observ- 
able A = ^ ^ ' input-independent in the perturbative regime at 

small X. gj^Q^ appears to be more input-dependent in the same range. The 
GRV-type parametrisations verify these properties. Other models, namely the 
BFKL kernel(QCD dipoles), DGLAP (hard pomeron singularity) give differ- 
ent predictions for A. At moderate Q^, there is a possibility of distinguishing 
these different perturbative QCD predictions in the near future. 
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I - The recently published 1994 results from HERA experiments on the proton 
structure function F2 have reached a high level of precision. It covers an extended 
kinematical range. In particular, it reaches very low values of x (x ~ 10~^) and 
(Q^ ~ 1.5 GeV"^). These data confirm with high statistics the strong rise of F2 when 
X becomes very small, first noticed in 1992 experiments [|I|. It is also observed in 
the newly reached kinematical range at small Q^. This rise has been quantified 
by a study of the observable 

^=(li!!i^)' 

where the brackets mean that A is obtained from a fit of the form F2 ~ at fixed 
and small values of a; (a; < .1). This observable has been previously proposed |^ 
as an interesting tool for discussing the various perturbative QCD expansions for 
the rise of F2. 

Two kinds of perturbative QCD predictions are available to explain this rise. 
At small X, the BFKL dynamics 0] naturally applies. It corresponds to the multi- 
Regge regime of perturbative QCD. In this approach, one sums up contributions of 
the type (^a^ln^j . In the present paper, we implement the BFKL dynamics using 
the recently developped QCD dipole model 0. It is based on the calculation of 
the infinite momentum wavefunction for arbitrary numbers of soft gluons in a heavy 
quark-antiquark (onium) state. Combined with fcy-factorization, this framework has 
been successfully applied to proton structure functions H, 0]. 



In another approach, the well-known DGLAP evolution equations ||I0| lead to 
alternative explanations for the rise of F2, based on the renormalization group evo- 
lution and the operator product expansion. These equations can then be expressed 
in term of moments in the cu-Mellin space. In this approach, the matrix elements 
of the local operators can be either regular at the right of = or singular in the 
cu-Mellin space. The first class of models was illustrated in the paper |jll| where a 
non-Regge behaviour for structure functions at small x was first suggested. In a sim- 
ilar framework, a perturbative evolution of valence-like input distributions lead to 
the parametrization of Gliick, Reya, Vogt (GRV) of structure functions [|T^, which 
gave satisfactory predictions for HERA. The second possibility, initiated some years 
ago and recently revived by Lopez, Barreiro, Yndurain (LBY), uses a singular input 
and also provides a satisfactory description of F2|T3|. 

Our aim is to study the properties of A (formula (0)) namely its dependence on 
the perturbative QCD origin of the rise of F2. The main results of this analysis are: 

i) In the GRV type parametrization, A is independent of the non-perturbative input 
for X < 5.10"^ and uniquely determined by the DGLAP kernel singularity. 

ii) The different types of perturbative QCD predictions, BFKL (dipole), DGLAP 
(GRV), DGLAP (LBY), are compared to the data on A. They are compatible with 
them but lead to significant differences at moderate (1 < < lOGeV^). This 
motivates precise measurements of A in this region. 
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II - Let us determine the predictions for A. First, we consider the direct and inverse 
Mellin transforms of the singlet structure function xxF^, namely 

Fsitu,Q^,fi^) = r dxx''-'xFs{x,Q\fx^) (2) 
Jo 

xF,{x, Q^ = / x^'^Fsico, Q^, fj.^)dcu, (3) 

where the integration line Reu = ujq is at the right of all singularities of xFs{uj, Q^, fi"^). 
In the DGLAP scheme xFs{u,Q'^, fi'^) and Fg{uj,Q'^, fi"^) (the gluon distribution in 
Mellin-moment space) verify 

where K{uj, Q^, /i^) is the Mellin transform of the DGLAP matrix kernel at the scale 
Q'^ ; The rightmost singularity of this kernel lies at uq = 0. Also is the 

Mellin transform of the singlet input which is either regular or singular at the right 
of lvq = 0, Fg{oj^^'^) is the Mellin transform of the gluon input. We then get two 
classes of models satisfying a DGLAP evolution ||14|| : 

(i) Either the rightmost i^-plane singularity is fixed by the DGLAP kernel singularity 
at = 0, we obtain the GRV type parametrization. 

(ii) Or we have the LBY type parametrization where the rightmost singularity lies 
at the right of = due to a singular input Fs{lj , Q"^ , fi^) , moreover its location is 
essentially not modified by the perturbative evolution. 

Let us observe that the BFKL dynamics also leads to a rightmost singularity at u 
greater than uj = 0. We will come back to this later. 

a -In the first case, we shall prove that the large behaviour of ^g^^f" is input 

independent and thus depends only on the DGLAP kernel. Let us start with a 
valence-like input, that is for which the a; = moment is well defined. In this case, 
the dominant a;— plane singularity is generated by the DGLAP evolution equations 

m 

Fs{uj, /i^) = f — ^exp z/+^ + — ^exp z/_^ ] (q + q) [uj) 



+ (exp z/+^ - exp z/_0 9i^) 

V+ — 



(5) 



where 



g, g, g are the valence-like input at the low scale = yU^. z/+, z/_, $^ are the first- 
order DGLAP kernels [|T^. A is the one-loop QCD scale and Nj is the number of 
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active flavours. The assymptotic form of i^+{uj) near tlie u; = singularity is 



UJ 



— a 



(7) 



where a is a constant (a = ^^|^ + |^ ~ 110/9 for 3 active flavors). Assymptotically, 



at a given Q near the uj = singularity we have 



Q , /X ) ~ ujfiuj) exp ( a)^ 



UJ 



(8) 



where /(cu) is an input- dependent function regular at = by hypothesis. Hence 
this fonction can be expanded as 



/M = /(0)[l + S,^*6,] 
Then, we can find an expression for xFs{x, Q^, fi"^) 



(9) 



1 rijjQ-\-ioc 

xFsix, Q\ /i') = — / u;/(o^)e "''^^/"e^^du; 

Z%T[ J LUn—ioO 



(10) 



I uiQ—ioo 

We now want to demonstrate that ^g'j^^fi is independent of the valence-like in- 
put. First, we can obtain an exact determination of this derivative when using the 
following property of Bessel functions 



1 
2i7r 



uo e 



n + l 



^ In 1 /x ' 



h[2{AN,i\nl/xY 



111 



Inserting the form (9) into (10), it is possible to compute xFs{x,Q^, fi"^) as was 
already done in [TH] for the lowest order in u. Defining lj and v 



UJ 



' In 1 /x 
V = 2{AN^^lnl/x)^, 



it is straightforward to see that 



d InxF, 



aini 



s - ^3 
UJ 



hiv) 



l + SA^-^(t;) 



(12) 
(13) 

(14) 



At this stage we can see that the only dependence on the input comes from the dif- 
ferent 6j, which can be shown to bring negligeable contributions. In our kinematical 
range for this study ~ 0.15 and x ~ 10~^, that is f ~ 10, uj ~ 0.4, we can use the 
assymptotic expansion for the Bessel functions (at large v) 



exp(t>j 
27rf 2 



1 - 



1 r(^ + i) 

2vT{v - I) 



(15) 
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which gives 



~ u- 



aini Uv) 



V V 



(16) 



So the absolute correction due to the term bi is of the order of which is neg- 
hgeable in the considered kinematical range. The terms bi,i > 1 are even more 
neghgeable. We deduce that ^^'j'^'fi does not depend on the different b'^ s. Hence the 

GRV type parametrizations ought to verify the following prediction 

(9 In- /2(f) 

which shows that this observable does not depend on the valence like input distri- 
bution. Using the precedent assymptotic expansion of the Bessel functions we can 
find an expansion (at large v) 

^'("' - 1 - A (18) 



which leads to 



Hence 



dlnxF, ^ ^ 5 

aini ~^\nl/x' ~ 41nl/x 

d"^ / 5 In xF. \ InxF, 



(19) 

This gives us an information on the concavity of the function X{Q^) in the DGLAP 
scheme with the rightmost singularity imposed by the DGLAP kernel. 



uses an 



b -In contrast to GRV type parametrizations, the LEY parametrization 
input c<j— plane singularity fixed in and located at the right of u = 0. Thus, our 
previous derivation does not apply in this case. Starting from the LEY formulation 
of the singlet structure function 

<e'> xF,{x,Q') = [Bs{Q')x-'' + Cs{Q')]{l - x)^(«') (21) 

where < > is the average charge for Nf flavours and Bs, Cs and u are Q^- 
dependent functions fl^. > defines the location of the rightmost cj— plane 
singularity and is Q^— independent, but for charm and bottom thresholds. At small 

^.x, (22) 

X 

up to the correction due to the phenomenological factor Cs{Q^). 

c -In the framework of the QCD dipole model, the EFKL dynamics also provides 

predictions for ^gl^f' ■ It gives ^ 

X 

Fs = Ca'/'x-"-^e-"^''''^ (23) 
Qo 



4 



where 

AaNc\n2 , /aNc 
ap = l + and a={ ^7(3 1x1- 24 

vr \ TT X J 

ap is the well-known BFKL Pomeron intercept which is a constant, since the strong 
coupling constant a is held fixed in this scheme. C and Qq are non-perturbative 
parameters to be determined by the fit. We get 

aini ^lni + 14^C(3)ln^i ^ ^ 

At fixed and x ^ 0, one recovers the usual BFKL Pomeron intercept. Note that 
a Q^-dependence stems from the last term in equation (pSj). Here we have, contrary 
to the GRV scheme 

9h^ > ° ^''^ 



III - Let us now discuss the phenomenological consequences of the previous cal- 
culations. In the following discussion, we have to take into account the differ- 
ence, due to the non-singlet contributions, between xFg and < F2, where 
< >= 2/9(5/18) for Nf = 3(4). However the non-singlet contribution is expected 
to be regular at = 0. Its QCD evolution receives no contribution from the gluon. 
Thus, one can apply our study of to F2 as well. In practice, our phenomenolog- 
ical discussion does include the non-singlet contribution when discussing the GRV 
and LEY parametrizations. Anyway, we have verified that this component is rather 
weak at small-x. 

At this stage, a comment is in order about the non-perturbative value of /x^. Indeed, 
in the leading-order version of the GRV parametrization, the evolution variable i^(/x^) 
(see formula (|])) is defined using the parameter values Nf = 3, fi = 480 MeV and 
^QCD = 232 MeV . We follow here an argument similar to the one given in ref. |Tl| . 



Starting from a perturbative scale Ql and evolving the structure function up to Q^, 
the renormalization group (DGLAP) predicts 

xFs oc e^«(«o)(^/(c^)) (27) 

where ^{Ql) is defined as in equation (|^) 

Note that ^ is the singular part of the DGLAP anomalous dimension i/+ at = 0. 

The form of the pref actor K~ ensures that the physical result does not depend on 
the arbitrariness of the factorization scale Qq in the perturbative range. The constant 
K cannot be predicted but has to be larger than one in order to get positive structure 
functions, using the functional form (|^), one can define /i^ in such a way that 

e(/i')=e(Qo)+lni^ (29) 
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where ^ is defined in equation (pHI). Since K > 1, one expect /i < Qo- Eventually, 
will be in the non-perturbative domain. This is indeed the case in the GRV 
parametrization, which uses an effective non-perturbative parameter /x^ and repro- 
duces the data. In conclusion, fi is an effective scale for the renormalization group 
evolution in the whole perturbative range Q » n without requiring the validity of 
the renormalization group in the non-perturbative range Q ^ fi. Note that a simi- 
lar argument also holds when applying the QCD dipole model for proton structure 
functions 0. 

In fig.l we compare the results of the GRV parametrization for ^^j^ with the 
prediction 

^h^ " 



This comparison is displayed in fig.l fonction of x for four different values 

of (5,20, 100,800 GeV'^). The agreement with the prediction ( ^0|) is reasonable 
for 10""^ < X < 2.10^^. Note that since this observable is a smooth function of x, 
the slope is extracted from a global fit of each bin by a function of the type 
We will come back to this point later. In fig. lb, we have plotted the Q^- 
dependence of ^g'^^F for x = 10~^ and 10~^. In both figures, a slight difference is 

observed between the two determinations. This might be due to the fact that the 
GRV parametrization is not an exact solution of the moment equations. A better 
agreement would be obtained with a slight change of /i^ from .23 to .24GeV'^. 
Fig. la and fig. lb illustrate our claim that the evolution of the considered observable 
is indeed dominated by the behaviour of the leading anomalous dimension. 
In fig. 2, we display the data on A defined in formula (|I|) and the prediction of the 
different parametrizations and evolution equations. The procedure for comparing 
models to data has been to use the same averaging in both cases. For each value of 
Q^, we have used the value of A given by a fit of the form x"'*' for 10~^ < a; < 10~^. 



Note that we have performed our own fit for E665 data [T^ while we used the 
published HI data for A with x < 10^^. We have checked that restricting the fit 
to the range 10~^ ^ x < 10^^ does not change the result. Fig. 2 shows that the A 
value obtained from the GRV parametrization is well reproduced by our prediction. 
It thus exhibits the universality property of A at small x in the whole range, 
since this universal value A is determined with essentially only one parameter /i^. 
Hence, a more accurate measurement of this observable may disentangle the nagging 
problem of the existence of the singular nature of the input. 

The different parametrizations are compatible with the data. The DGLAP (LEY) 
parametrization is displayed for > 10 GeV^ in its domain of validity. We note 
in this case that A is slowly varying with and lower than the A^ = .36 in formula 
(p2D . This might be due to the extra singlet factor Cs in equation (^1]) . 
The BFKL (dipole) prediction is satisfactory. It is interesting to note that it predicts 
a variation of A, showing that a model satisfying the BFKL dynamics combined 
with fc^-factorization induces a variation of the slope, even if the tu— singularity 
is fixed. We note a difference with the DGLAP predictions in the small range 



6 



{I < Q"^ < 10 GeV"^). Accurate data in this range might distinguish between the 
two different approaches. 

In this context, we have also included in fig. 2 the very recent HI 1995 preliminary 
results for A [0. They unravel a tendency in favour of the DGLAP (kernel) scenario. 



This new measurement, if confirmed and made more precise in the future, shows 
the reliability of measuring A for distinguishing the QCD origin of the rise of F2. 
The quality of the experimental determination deserves more study at the theoretical 
level, for instance a more precise determination of the QCD dipole predictions at low 
Q^, a study of the flavour threshold effects in the LEY approach and a DGLAP kernel 
considered at higher perturbation order. Work is in progress in those directions. 

One could wonder whether the observable ^gi^^Q2 (or^i^|^) obeys the same uni- 
versality property as in GRV type models, since the DGLAP kernel is the 

dominant singularity in the cu— plane. The following study shows that this is not 
the case. It turns out that it depends strongly on the input, and it is thus not 
dominated by the anomalous dimension. A developpement similar to the previous 
one for the determination of ^J^^^" leads to 

aml/x 

dlnxF, ^ AN,h{v) + ^,hu'h+^{v) 
with v+ioj) ~ — . Now, if we take an accurate developpement for z/+, see (7), we 



get 



The assymptotic expansion for the Bessel functions leads to 

Note that the main contribution to ^ ^^m^" behaves as 4 in contradistinction with the 

previous case where gj^"^^ is of order uj. This is the origin of the non universality 
of this observable as shown below. 

Indeed the term hi contributes to with a coefficient to be compared with 

its contribution to IHr^ where the coefficient is — . As - 0.4, v ~ 10, AN^ = 12, 
the correction in bi to the kernel contribution 

a+ _ 34 

d4 u l2{v) 

is important, so this observable is very sensitive to the corrections due to the terms 
bi- Moreover ~ a, so this derivative is sensitive to an higher order of de- 

veloppement of i^+{u!) in u. We can undestand this point by writing exactly the 
derivative with respect to C, 

dlnxFs 1 f^o+ico ^ 

r,, = 77- / ujf{u)u4u)e "'"^/"e - ^du, (35) 
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The terms coming from the development of i^+iuj) will be mixed with those from 
f{uj) and we can not separate both contributions as it is the case for . Thus 

^^gl''^" does not obey the same universality property as §{^^^- We can notice here 
that another explanation using a saddle point method introduced in ref . |^ leads by 
a different method to the same conclusion [Q. 

To summarize, we have shown that, provided that the non-perturbative input is 
regular at the right of the u = singularity of the dominant DGLAP anomalous 

dimension, the observable A = ^ is input-independent in the perturbative 

q2 j.ggij2ie at small x. Other models, namely BFKL(dipole), DGLAP (LBY) give 
different values for A, which are compatible with present published data. This is an 
incentive for the experimentalists to get a better accuracy, and for the theoreticians 
to refine the predictions in order to distinguish these different QCD interpretations 
in the near future. 
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FIGURE CAPTIONS 

Figure la Comparison of tlie GRV parametrization for ^g^^f" witli our predic- 

X 

tion ^-^^ for the same observable (called DGLAP (Kernel) on ttie figure). This com- 
parison is displayed as a fonction of x for 4 different values of (5, 20, 100, 800 GeV^) 

Figure lb Comparison of the GRV parametrization for ^g^^f" with our prediction 

cuj^l^ for the same observable (called DGLAP (Kernel) on the figure). This com- 
parison is displayed as a fonction of for two different values of x laying in the 
kinematical range of our study : x = 10~^,a; = 10"''. 

Figure 2 Display of the data on A (HI: 0, |l^ and E665: |jl6|) compared with the 



prediction of the different parametrizations and evolution equations. 
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